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by 

D. P. Gaver 

v 

_and 

-V 

P. A. Jacobs 

Operations Research Department 
Naval Postgraduate School 
Monterey, CA 93940 

1 . Introduction and Assumptions 

The integrity of a great many physical structures is 
potentially threatened by combinations of physical loads of vary- 
ing magnitudes from various sources. We are thinking of structures 
such as buildings (for instance those that house or contain 
nuclear power plant elements) ; aircraft and spacecraft; electrical 
transmission networks; bridges, piers, and dams; and offshore 
oil drilling rigs that experience loads from wind, snow and ice, 
tides, earthquakes, and so forth. In many instances the total 
load or stress experienced by a structure varies in time in an 
apparently random fashion. Certain load components vary rather 
slowly; for example, that component resulting from snow and ice 
accumulation; others occur more nearly as impulses, such 
as those associated with winds or earthquakes. The problem is 
to design structures to withstand the superposition of loads 
from many sources with at least an approximately understood 
(high) probability. In engineering terms we wish to work towards 



1 



developing a rational safety factor technique for designing 
structures to withstand the combination of loads anticipated. 

The purpose of this paper is to describe and investigate certain 
simple but somewhat realistic probabilistic load models for 
use in design, and perhaps safety, assessment of structures. 

In this paper we confine attention to the superposition 
of just two load types: shock loads , and constant loads . For 

example, wind gusts, flash floods, and earthquakes have varying 
magnitudes and have relatively short durations in comparison to 
the times between their occurrences; these will be modeled as 
instantaneous shock loads. On the other hand snow, ice, or 
water accumulation, or even the presence of slow-moving vehicles 
or furniture, present loads that remain nearly constant in time, 
occasionally changing to new levels; these will be modeled as 
constant loads that change infrequently. Throughout this investi- 
gation it will be assumed that the effective stress exerted by 
several types of loads acting simultaneously can be expressed 
as a linear combination (actually, sum) of those component loads, 
the loads being treated as stochastic processes. Somewhat 
similar load combination models have been studied in the past; 
see the work of Bosshardt [1975], Larrabee and Cornell [1978], 
McGuire and Cornell [1974], Pier and Cornell [1973], Weh [1977], 
and others. Note that the present study utilizes the notion of 
stochastic point processes as a modeling tool; see Cox and Miller 
[1965] for an introduction, and Lewis (ed.) [1972] for a more 

extensive treatment. 
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1.1. Modeling Assumptions and Problems 



Assume that the successive magnitudes of shock loads are 
independently and identically distributed (i.i.d.) random variables 
with common continuous distribution function G(y); G(y) = 0 
for y _< 0 . The shock loads' appearance is regulated by a Poisson 
process with rate y. Furthermore, the constant loads have i.i.d. 
magnitudes with distribution F(x); F(x) = 0 for x ^ 0. Constant 
loads change at moments of a Poisson process with rate A . 

Finally, the shock and constant load processes are, for the 
most part, taken to be statistically independent, although 
this assumption can be relaxed at times without difficulty. 



Let X(t) (respectively Y(t)) be the magnitude of the 
constant (shock) load process at time t; (Y(t) will usually 
be zero) . Put Z(t) = X(t) + Y(t), the superposition of the two 
loads at time t, and M(t) = sup g < t Z(s) , the maximum load 



combination in [0,t]. See Figure 1. Let T y ={inf t_>0:Z(t) > x) , 

CONSTANT LOAO PROCESS 

x ( t ) 1 



Y(!) 



SHOCK LOAD PROCESS 



Z ( t ) 



LOAD COMBINATION PROCESS 



Fig. 1 

THE LOAD COMBINATION 
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i.e. the first-passge time to a combined stress level x. Unless fatigue 
or some other cumulative effect is in operation T will represent 
the time to failure of a structure whose strength is x, and that 
is subjected to a stress history (Z(t) , t >_ 0) . 

Section 2 is concerned with the description of the distri- 
bution of the maximum process and the distribution of T . The 
first step is to obtain Laplace transforms with respect to time. 
Asymptotic results (x + °°) for E [T ] and the distribution of 

X 

the normalized value T* = T / (E [T ]) will then be presented. In 

particular it will be shown that the distribution of T* is 

approximately a unit exponential distribution as x becomes large 

and so does E [T ] . 

x 

Section 3 is devoted to a study of the asymptotic 
properties of the "maximum process," M(t), as t °°. These 
results relate to those of Welsch [1972] and O'Brien [ 1974a, b]; 
they are seen to extend the classical "extreme value" results 
of Gumbel [1958] and Gnedenko [ 1943]. we present results for 
the asymptotic distribution of M(t) and for the joint distri- 
bution of the first and second maximum load combination in 
[0 , t ] as t -*• °°. 
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2 . 



The Laplace Transform of the Distribution of M(t) 



and a First-Passage Time Limit Theorem 



2.1. Towards the Distribution of M(t) 



Consider first the distribution function of the maximum 



process, H (t) = P{M(t) £ x} , for x > 0. Its Laplace transform 
is available immediately by noting that if T^ is the time of 
the first change in the magnitude of the constant load process, 
then 

H (t) = P{M(t) £ x, T x > t} + P{M(t) £ x, £ t} 



e ^ / exp{-yt G(x-y)} F(dy) 



0 



l— _ "l 

+ / Ae V dv / exp{-yv G(x-y)} H (t-v) F(dy) . (2.1) 

« « X 



o o 



Next take Laplace transforms with respect to t: 




H (t)dt = M (£) + AM (?) h(C) ; 

X X XX 



( 2 . 2 ) 



reversal of the order of integration provides that 



x 



-1 



M ( O = } U + A + yG(x-y) ] 
x 0 



F(dy) . 



(2.3) 
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Hence, 



M U) 

V 5) - r --\ru ■ 

It seems to be difficult to solve (2.1) (invert (2.4)) in any 
simple convenient form for any interesting choice of the distri- 
butions F and G. However, useful information can still be 
gleaned from (2.4). First, it is clear that if T is the 
first-passage time to combined load x, so 



T = inf(t > 0:M(t) > x} , 

x — 

then, since P{M(t) £ x} = PtT^ > t} , 

oo 

h (?) = / e"^ fc P { T > t} dt 

x o x 

Let £ -*■ 0 in (2.4) to find that 



m(x) = E [T ] = h (0) 
x x 



M (0) 

r - AM (0) 



(2.5) 



( 2 . 6 ) 



(2.7) 



We now record some expressions for the mean first-passage time 
to x when specific distributions for shock and constant load 
maqnitudes are in force . 
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— x 

Example 2.1 . Identical exponentials: F(x) = G(x) = 1 - e , x _> 0 . 

In this case 



... 1 -x, yxe ye „ A + ye 

M x (0) = X (1 - e ) - — ^ — In [xmr - J 



and 



E [T 



A(l-e x ) - yxe x + ye x £n [ ( A + y ),/ ( A + ye x ) ] 



1 + x + £n [ ( A + y) / ( A + ye x ) ] 



(2.3) 



- (-) 
y \ x / 



as x -+ 



Example 2.2 . Different exponentials: 

F ( x) = e -ax , G ( x) = e kx where ^ = 2. 

In this case 

M (0) = / [A + ye -k(x-y) ] _1 ae _ay dy 
x 0 

= A' 2 e -2kx N (x) 

where 

N(x) = A [-1 + e 2kx ] + 2y[l-e kx ] 

9 -ley -1 

+2y [-kx + £n [ ( A + y ) ( A + ye ) ]]. 
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Thus 



Et V 



A 1 N (x) 



[Ae 2kx - N (x) ] 




(2.9) 



2 

Example 2.3 . F(x) = e x , G(x) = e x . In this case 



M x (0) = j [1 - e -x ] + ^ e" X N ( x) 

where 

N (x) = (e _X - 1) [(A + y) (A + ye~ X ) ] _1 

Thus 



E [T ] 
x 



Y e X [1 - e X + ye x N (x) ] [1 - yN (x) ] ^ 



A + y x 
A(A+y) + y e 



Example 2.4 . 



( 2 . 10 ) 



F (x) = 



a l X 



-a 



1/a 



x < a 



x 1 a i /a / 



and 



G (x) = 



a_x 
z 



-a 



x < a 



x > a 



1/a 

2 

1/a 



a l' a 2 



where a > 0 and 



> 0. In this case 



as x -»■ oo 



Example 2.5 . 



F ( x) = 



and 



G (x) = 



-a 



x 



-3 



if x < 1 , 

if x ^ 1 , 

if x < 1 , 

if x > 1 . 



If a < 3 , then as x 



E [T ] ~ y x a 

X A 



If 3 < a, then as x 



[T ] ~ - x 6 
x J y 



( 2 . 12 ) 



The next result is a limiting result for the first 
time the load combination process exceeds a given level x. 



THEOREM (2.1). The limiting distribution of T^ is exponential, 
in the sense that 



lim P{m(x) ^T > t) = e 
x 

x-x 0 

where 

x Q = inf (t:F*G(t) = 1} . 

and 

m ( x) = E [T x J . 
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Proof. Put p = pA ^ and rearrange (2.3) so that the 
denominator of (2.7) is in the form 

x _i 

1 - AM (0) = F (x) + / pG(x-y) [1 + pG(x-y)] F(dy) = h^x). (2. 

x 0 

Now apply the bounded convergence theorem and the fact that F 
and G have densities to show that m(x) 00 as x -»• x^ . It 
next follows by a change of variables in (2.4) that 



/ P(m(x) ~ T > t}dt = m(x) h (y(x)) 

0 xx 



M x (y (x) ) 

m (x) [1-M (y (x) ) ] ( 2 

X 



where y(x) = £m(x) 
to (2.3), 



-1 



By the bounded convergence theorem 



lim M x (y (x) ) = - . 

x>x 0 



Again from (2.3) we have 



applied 

(2.15) 



1 - AM^ (y (x) ) 
x 

= F ( x) + / [y (x) A 1 + pG(x-y) ] [y (x) A -1 + 1 + pG(x-y)] -1 F(dy) . 

0 



Hence , 



m(x) [1 - AM (y (x) ) ] 

X 

= Y / [y(x)A 1 + 1 + pG(x-y) ] 1 F (dy) 

A 0 

x _ 

+ m(x)[F(x)+ / pG(x-y) [y(x)A + 1 + pG(x-y) ] F(dy)] . 

0 

Let 

— ^ _ _ I 

h~(x) = F ( x) + / pG (x-y) [y (x) A + 1 + pG(x-y)] F(dy) ; 

0 

then 

h-^ (x)- h 2 (x) 

x _ 

= / F (dy ) pG (x-y) {-[l + pG(x-y)]+ [y(x)A +l+pG(x-y)]} 

0 

*{[1 + pG (x-y) ] [y(x)A _1 + 1 + pG(x-y)]}" 1 . 

X _ 

= / F (dy) pG (x-y) y(x)A {[l+pG(x-y)] [y(x)A _1 + l+ C(x-y)]}" 1 . 

0 

= Y (x) A 1 k (x) 
where 

x _ _ _i _ 

k(x) = / pG (x-y)([l + pG (x-y) ] [y (x) A +l + pG(x-y)]} F(dy). 

0 
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By the bounded convergence theorem 



lim 
x x 



0 



k ( x) 



0 . 



Recall that 



Hence , 



m(x) 



j + 0 (x) 
h 1 (x) 



h^ ( x) - h 2 (x) 
m(x) 



£h^ ( x) A 
A ^ + 0(x) 



k(x) 

h^TxT 



K 

1 + 0(x) 



k (x) 



which tends to zero as x x Q . Therefore, 



lim m(x) [1 - AM (y(x) ) ] = £A ^ + A ^ lim 

x x o X x * x 0 



1 + 



h 2 (x) - h^ (x) 
h”( x) 



= U+ 1) • A 1 . 

The result now follows from (2.14), (2.15) and the unicity 

result for Laplace transforms. 
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3. 



Asymptotic Results for the Maximum Load Combination 



In this section we will study the asymptotic behavior of 
the maximum load combination to occur during time interval ( 0 , t ] 
as t Most results concerning maxima of non-Gaussian random 

variables are for variables in discrete time. Thus we will first 
obtain results for an embedded discrete time load combination 
process . 

Let be the time of the nth change in magnitude of 

the load combination process; that is. 



S 1 = inf { t > 0 : Z ( t) f Z(t-)} 



and for n > 1 

S n = inf {t > S n . 1 ;Z(t) ? Z ( t- ) } . 



(3.1) 



The change in the load combination at time S n may be due to 

either a change in the constant load process or to the arrival 

of a shock load. Let Z = Z(S ) (respectively, X = X(s ), 

n n * J n n 

an< 3 = Y(S n )) be the magnitude of the load combination 

(respectively constant load and shock load) at time S . Put 



M = max (Z ) 
0 < k < n K 



(3.2) 



Note that the times {S. n > are the arrival times of a Poisson 

process with rate A + y and are independent of { Z } . Further 

n ' 
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